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Abstract 



The phases in a generic low-energy supersymmetric model are severely constrained 
by the experimental upper bounds on the electric dipole moments of the electron 



and the neutron. Coupled with the requirement of radiative electroweak symmetry 



| breaking, this results in a large degree of fine tuning of the phase parameters at the 

0> unification scale. In supergravity type models, this corresponds to very highly tuned 

■ values for the phases of the bilinear Higgs coupling parameter B and the universal 

r S^ • trilinear coupling Aq. We identify a cancellation/enhancement mechanism associated 
Oh' 

with the renormalization group evolution of B, which, in turn, reduces such fine-tuning 

Oh! 

D , quite appreciably without taking recourse to very large masses for the supersymmetric 



partners. We find a significant amount of reduction of this fine-tuning in non-universal 
gaugino mass models that do not introduce any new phases. 



*U ■ PACS numbers:13.40.Em,04.65.+e,12.60Jv,14.20.Dh 



1 Introduction 

Low energy supersymmetry (SUSY) [1] has been playing a central role in the quest for 
physics beyond the standard model (SM). Since phenomenological consistency requires SUSY 
to be broken, and broken softly (so as not to reintroduce any quadratic divergence), the 
Lagrangian of the Minimal Supersymmetric Standard Model (MSSM) [2,3] includes soft and 
gauge invariant SUSY breaking terms. While the generic MSSM Lagrangian may contain 
many arbitrary soft terms, specific models for SUSY breaking have been proposed that 

1 Emails: tpuc@iacs.res.in, debchou@physics.du.ac.in, tpdd@iacs.res.in 
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provide relationships between the MSSM parameters. Incorporating well-motivated new 
interactions and particles at high mass scales, such scenarios drastically reduce the large 
number of unknown parameters in the MSSM to only a few, thereby making the model more 
predictive. We will focus here only on supergravity (SUGRA) [4, 5] type of models where 
SUSY is considered as a local symmetry. These models incorporate a hidden sector wherein 
SUSY is broken, and a visible sector where the MSSM fields reside and to which the breaking 
is communicated by gravitational interactions. In N = 1 SUGRA, which incorporates grand 
unification, one has a choice of three functions in building a model [4-6] , namely the gauge 
kinetic energy function f a p(zi), the Kahler potential K(zi, z\), and the superpotential W(zi), 
where Zi refer to matter fields. In mSUGRA, the minimal version of the model, one has a 
flat Kahler potential and a flat gauge kinetic energy function. The corresponding soft SUSY 
breaking sector is characterized by only a few parameters, normally specified at the scale 
of the grand unified theory (GUT) viz. Mq ~ 2 x 10 16 GeV [7,8]. These are the universal 
gaugino mass mi, the universal scalar mass m , the universal trilinear coupling A and the 

2 

universal bilinear coupling B . In addition to these, there is a superpotential parameter, 
namely the Higgs mixing term /z . Unlike in the SM, where the breaking of the electroweak 
symmetry necessitates the explicit introduction of a negative valued scalar mass-squared, 
in a generic SUGRA model, the said breaking can be realized even for a positive mass- 
squared term in the bare Lagrangian, thanks to radiative corrections [4]. In other words, 
the renormalization of the soft SUSY breaking terms as one moves from the unification scale 
down to the electroweak scale automatically engenders a negative mass-squared thereby 
breaking the symmetry [9-12]. In a similar vein, the low energy parameters of the MSSM 
(which are quite large in number) are obtained from only a few unification scale parameters 
via the renormalization group equations (RGE) [12] integrated from Mq to the electroweak 
scale (~ Mz). The two minimization conditions for the Higgs potential then eliminate \xq 
(except for its sign) on the one hand, and, on the other, relate Bq to tan/5 (= (Hu) / (Hp)), 
the ratio of Higgs vacuum expectation values. Thus mSUGRA may be characterized by 
tan/3, mi/2, mo, A and sign(/i) 2 . With all the low energy parameters of the MSSM being 
generated in terms of these few parameters, one has a considerable amount of predictivity 
for the MSSM spectrum. 

A different problem remains though, namely that of the SUSY CP violating phases. Many 

2 Our choice of sign for \i and Aq follows the standard convention of Ref [13]. 
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phases of SUGRA models can be rotated away. In an universal scenario like mSUGRA, the 
gaugino masses can be considered real with the result that only two combinations of phases 
(beyond the Cabibbo-Kobayashi-Maskawa quark mixing (CKM) phase already present in 
the SM) are physical. A convenient choice for the two is given by 4>a f° r (at Mq) and 9 b 
for the £>-parameter at the electroweak scale. It should be noted though that many analyses 
prefer to work with 6^, the phase of n, instead of 9b- An advantage of this latter choice is 
that 9^ Q ~ 6^ since 9^ does not run up to the one-loop level. These different descriptions 
can be understood in terms of U(1)r and U(1)pq (Peccei-Quinn) symmetries and the choice 
of reparametrization invariant combinations of phases, a discussion of which may be found 
in Refs. [3, 14]. A selection of past analyses using 9b as an input parameter may be seen in 
Refs. [15-20]. Here we note that a choice of 9b instead of 9^ as a phase parameter makes 
the entire set of input parameters to be of soft-breaking origin. 

A few important points need to be noted in the context of the SUSY CP problem. 
The latter arises from the fact that the phases are highly constrained by the experimental 
limits on the electric dipole moments (EDM) of the electron and the neutron [14-19,21,22]. 
Consequently, we are forced to admit one of the three eventualities: 

1. The phase 9b is very small — O(10 -2 ) or 0(1CT 3 ) — if the superpartners are not consid- 
ered to be very heavy 3 . In addition, the phases of the A-parameters at the electroweak 
scale are also constrained. In mSUGRA with phases, the requirement of having a 
very small 9b typically translates into a relatively large but highly fine-tuned value for 
arg(_B ) (i.e., B at Mq). This, in turn, constrains the phase 0a o of A , although to a 
somewhat lesser degree. The fact that the issue of fine-tuning in phases at the GUT 
scale arises out of the combined requirement of satisfying the EDM constraints and the 
radiative electroweak symmetry breaking was discussed in great detail in Refs. [15-17] 
as well as in Refs. [18, 19]. In this paper we try to focus our attention on this problem 
by looking at suitable models beyond mSUGRA that can have unique features in the 
evolution of B. 

2. The phases are large and less fine-tuned but the sparticles are massive. Of course, 
fully ameliorating the SUSY CP problem in this fashion requires that the sfermions 
be super-massive, thereby aggravating the problem of the little mass hierarchy in the 

3 9b may reach up to ~ 0.1 in the focus point zone [23]. 
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Higgs sector. We will investigate whether the amount of fine-tuning can be reduced 
even while one considers a lighter sparticle spectra. 

3. Finally there is the possibility that the SUSY breaking parameters may have special 
pockets where there can be a large amount of internal cancellations between the dia- 
grams contributing to the electric dipole moments of electron and neutron [21]. This 
means that phases could be large while sparticle masses are significantly light. This 
scenario is highly parameter dependent and clearly depends on very delicate cancel- 
lations. Hence we will not include this in our work while trying to focus on generic 
behaviors. 

As mentioned above, we would like to address the first and the second issues in this 
analysis. We are particularly interested in exploring the possible role of non-universal gaugino 
masses (NUGM) in reducing the fine-tuning in the phase 9b - To quantify the latter, we 
consider a naturalness like measure of the form 

$ = [A9 B JA9 B ]e B ^. (1) 

A large value for $ would mean a lesser degree of fine-tuning of 6b with respect to a variation 
in 9b satisfying the EDM constraints. The phase-derivative is evaluated at 9b ~ with the 
choice being dictated by the fact that the EDM constraints force \9 B \ to be close to zero. 
Thus, this is a restrictive definition compared to the type of fine-tuning defined in Ref. [16]. 

We will see that the issue of such fine-tuning of phase can be addressed by focusing on 
scenarios where there is a large evolution of the bi-linear Higgs coupling parameter B between 
the electroweak scale and the GUT scale. The evolution of B depends on the U(l) and the 
SU (2) gaugino masses, the trilinear couplings and tan j3. Within mSUGRA, in addition to 
the evolution of \B\ being typically small, the phase 9b also turns out to be quite fine-tuned 
(i.e. <3> tends to be small). In other words, for a given 9b satisfying the EDM constraints, the 
variation A9b that still is consistent with the constraints is generally much smaller than the 
variation A9b allowed at the electroweak scale [15]. As we will see, the evolution in \B\ may 
be enhanced by appropriate mass relationships between the gauginos that are away from 
universality at Mq. At the same time, these would help in reducing the above-mentioned 
fine-tuning so that $ can be significantly increased in specific NUGM scenarios. 

We, however, desist from choosing an arbitrary non- universal gaugino mass scenario since 
that will introduce new phases [17]. As we will see in Sec. 2, non-universalities in gaugino 
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masses may originate from a non-trivial gauge kinetic energy function. The latter is a 
function of chiral superfields and transforms as a symmetric product of the adjoint represen- 
tations of the underlying gauge group. This leaves f a p with the possibility of being in one or 
more of several representations, one of which is the singlet. While the choice of the singlet 
corresponds to mSUGRA, the non-singlet representations give rise to non-universalities in 
the gaugino masses. It is possible to identify a suitable non-singlet representation in isolation 
(i.e., we will not combine a non-singlet representation with the singlet or other non-singlet 
representations) whose gaugino mass pattern is effective in generating a large evolution in 
B. At the same time, there will be no additional phases to worry about since the overall 
phase of the gaugino masses can be rotated away in a fashion similar to that in mSUGRA. 

In this paper, we will analyze the consequences of a large evolution of the 5-parameter 
(mostly in the presence of such non-universalities) on the CP violating phases. Here, the 
basic input parameters are tan/3, mo, mi (providing with definite NUGM patterns), \Aq\ 
along with its phase 4>a and the phase 9 b of B given at the electroweak scale (~ Mz)- Note 
that \B\ at the electroweak scale is obtained via radiative electroweak symmetry breaking 
(REWSB) condition. Subsequently, \B \, the GUT scale magnitude of the S-parameter along 
with its phase 9b is obtained via RGEs. We will identify broad but correlated regions of 
parameter space where there can be a significant degree of reduction of the phase sensitivity 
while going from mSUGRA to a type of NUGM models. 

The paper is organized as follows. In Sec. 2, we discuss the non-universal gaugino mass 
models. The study of the relevant contributions from different sectors in the associated RGEs 
of B and A parameters allows us to identify the non-singlet representations which provide 
with a large evolution in B. We will probe the parameter space that is suitable for reducing 
the amount of fine-tuning in the CP violating phases. In Sec. 3, we present the numerical 
results for the evolution of B. An analysis in the absence of phases points us to the favored 
regions of parameter spaces. On inclusion of phases, this facilitates the identification of the 
regions with significantly reduced level of fine-tuning, Finally, we conclude in Sec. 4. 
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2 Non-universal gaugino masses and enhanced evolu- 
tion of B 



Non-universality in gaugino masses may originate from a non-trivial gauge kinetic energy 
function f a p which, in turn, is a function of the chiral superfields in the theory. The indices 
a, (3 run over the generators of the gauge group (for example, a = 1,2, ...24 for SU(5)). 
The gaugino mass matrix is given by 

M aP = \e G ^G a (G- 1 ) b a (df: i /dz* b )f;, 1 (2) 

where G = — \n[K 6 WW*\ — k 2 K. Here, W is the superpotential, K(z,z*) is the Kahler 
potential, z a are the complex scalar fields, and k = (8ttGn)~^ = 0.41 x 10~ 18 GeV" 1 with 
Gn being Newton's constant. The functions f a p may have non-trivial field contents, or 
in other words, may contain combinations of field transforming as either singlet or non- 
singlet irreducible representations [24]. With the gauginos being Majorana particles, f a p, of 
necessity, must be contained in the symmetric product of the adjoint representations of the 
gauge group. For example, in the case of SU(5), 

f a p D (24 ® 24) aym = 1 © 24 © 75 © 200 . (3) 

For the singlet case, one has f a p = 5 a p which indeed leads to universality of gaugino masses. 
Similarly, the non-singlet representations will give rise to non-universal gaugino masses. 

In general Mi(Mc) = mi ^ r C r n£, where CV's give the relative weights of each contribut- 
ing representation and n[, for the subgroup i, are essentially the Clebsch-Gordan coefficients 
corresponding to the breaking by the adjoint Higgs field [24-26]. For the case of 577(5), 
the coefficients n\ are displayed in Table 1. Clearly, the non-singlet representations have 
characteristic mass relationships for the gaugino masses at the GUT scale. Past analyses 
exploring various phenomenological implications of such non-universality may be found in 
Refs. [24,25,27-29]. 

As we shall argue later, the adjoint representation r = 24 for f a @ (NUGM:24 in the 
notation of Table 1) is the most interesting one in the context of the present investigation. 
Consequently, we will analyze this case in isolation, or, in other words, assume that the sole 
contribution to f a p is from a 24-plet structure. Apart from reducing the number of free 
parameters, this has the additional advantage that no new phase degree of freedom for the 
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gaugino masses is introduced. With the gaugino mass ratios at the GUT scale now being 
given by M 3 (M G ) : M 2 (M G ) : M\(Mq) = 1 : -3/2 : -1/2, for a positive gluino mass, the 
other two gaugino mass parameters are negative, a signature different from mSUGRA. This 
indeed would turn out to be useful in our quest. As mentioned earlier, we only consider 
either C\ = 1 (mSUGRA) or C 2 4 = 1 (NUGM:24) with all other C r 's assumed to be zero. 



r 


Label 


M 3 G 


M 2 G 




1 


mSUGRA 


1 


1 


1 


24 


NUGM:24 


2 


-3 


-1 


75 


NUGM:75 


1 


3 


-5 


200 


NUGM:200 


1 


2 


10 



Table 1: The coefficients n\ as pertaining to the SU(3), SU(2) and U(l) gaugino masses at 
the GUT scale for different representations of SU(5). 

An analogous analysis with SO(10) as the underlying gauge group is also possible [29,30], 
though we will not investigate it in this paper. Similar to Eq.3 here, one has (45 x 45)^^ = 
1 + 54 + 210 + 770. If the symmetry breaking pattern is 5O(10) SU(4) x SU{2) x SU{2) -> 
SU(3) x SU{2) x U(l), one finds from the 54-plet that M 3 (M G ) : M 2 (M G ) : M^Mg) = 1 : 
—3/2 : —1. This pattern is quite similar to NUGM:24 as can be ascertained from Table 1. 
We would like to comment at this point that, in general, such non- universal gaugino mass 
scenarios change the gauge coupling unification conditions [24,26]. However, it is still possible 
to find specific conditions [24,31] under which the usual gauge coupling unification condition 
remains unaltered and we consider this in our work. Note though that our results are quite 
robust and have very little dependence on the exact details of the spectrum. 

2.1 Nature of evolution of B with real parameters 

We now identify the differences between mSUGRA and NUGM:24 in regard to the evolution 
of the 5-parameter in the absence of CP violating SUSY phases. This, in turn, will help 
us in understanding the evolution of 9b upon the inclusion of the phases (see Refs. [17-20] 
for past analyses discussing phase evolutions). Note that /i 2 and B are determined via the 
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REWSB condition, viz. 

I ,2 1 M 2 m H D ~ m H v ^n 2 (3 Ei - E 2 tan 2 

1/11 2 z tan 2 /3-l tan 2 /3-l ( 4 ) 

sin(2/3) = 2|5/x|/(m| D + m\ + 2/i 2 + Ei + E 2 ) , 

where Ej represent the one-loop corrections [32,33]. The Higgs scalar mass parameters rriH D 
and rriHjj, and thereby /i 2 and B depend quite strongly on m as well as on mi. To one- loop 
order, the running of the B parameter has two additive components, the first proportional 
to the gaugino masses and the second depending on a combination of the trilinear couplings 
and the Yukawa couplings [12,17], namely, 
dB 3 

— = (3a 2 m 2 + + (3Y t A t + 3Y b A b + Y T A T ) , (5) 

at 5 

where t = 1?i{Mq/Q 2 ) with Q being the renormalization scale. 6l\ = o:j/(47r) are the scaled 
gauge coupling constants (with ot\ = fay) and m; for i = 1,2,3 are the running gaugino 
masses. Furthermore, Yi represent the squared Yukawa couplings, e.g, Y t = yf/(47i) 2 where 
y t is the top Yukawa coupling. In a similar vein, the evolution of the trilinear terms is given 
by 

dA t ,16 13 , . „ A 

—— = -(—a 3 m 3 + 3a2m2 + —aim 1 )-6Y t A t -Y b A b 
at 3 15 

—jT = -(^-^3^3 + 3a 2 m 2 + — c^mi) - F t A t - 6Y b A b - Y T A T (6) 
at 3 15 

= -(3a 2 m 2 + ^aimi) - 3Y" 6 v4 fe - AY T A T . 
at 5 

For small tan j3, the contributions from the bottom quark and tau Yukawa couplings y b 
and y T may be neglected, and the RGEs approximately integrated to obtain [15] 

B-Bo^ 00 ®' 1 Ao-C(t)mi , (7) 
z 2 

where Do(t) = 1 — 6Y(t)F(t)/E(t) with t corresponding to the electroweak scale. The 
functions E(t) and F(t) encapsulate the running of the gauge coupling constants, viz, 

E(t) = (1 + /3 3 t) 16/(3b3) (l + (3 2 t) 3/b2 (l + At) 13/(15fel) 



F(t) = f E(t')dt' 
Jo 



where $ = 6^3:^(0) and {bi,b 2 ,b 3 ) = (33/5, 1, —3) are the coefficients in the respective one- 
loop beta-functions. Of course, unification imposes the boundary condition that «j(0) = 
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«g ~ 1/24. At the top mass scale (Q = m t ), D — 1 — (m t /200 sin/3) 2 ^ 0.2 is indeed a 
very good approximation. The function C(t), in Eq.7, on the other hand, is given by 

C(t) = -1(1 - D )^ + Uh 2 + ^ , (8) 



where 

Kit) 



t 



+ 

H 3 (t) = [ E(t')H 2 (t')dt' 
Jo 

H 2 (t) = a(0) (^-h 3 + 3h 2 + ^h 1 
\ 6 15 

For the generic (NUGM) case, the above results remain the same except that [25] 

K{t) — k(t) = hi(t) ^ . (9) 

mi 

2 

Note that, in dB/dt, the gaugino contribution is positive for mSUGRA, but negative for 
NUGM:24. Thus, it is useful to understand the nature of evolution of trilinear couplings in 
either scenario so as to evaluate their role in the evolution of B. For the mSUGRA case, the 
gaugino contributions to dAi/dt are always negative (vide Eq.6). Hence, it is obvious that if 
Aq not be too large, then Ai would typically turn negative by the electroweak scale. In fact, 
the large gluino contributions render both A t and Ab negative well above the electroweak 
scale. This implies, that in this case (mSUGRA), the two pieces in dB/dt would tend to 
cancel each other, an effect also manifested by the smallness of C in Eq.7. In turn, this leads 
to a small value for AB = \B$ — B\ in mSUGRA. 

Comparing the evolution of the trilinear terms in NUGM:24 with that in mSUGRA, 
it turns out that a qualitative difference arises only in the case of A T , while for A t and 
Ab the difference between the scenarios is only a quantitative one. This is easy to under- 
stand given the overwhelming dominance, in the last two cases, of the gluino contribution 
over those from the electroweak gauginos. Specifically, for Aq = 0, A T at the weak scale 
comes to be negative for mSUGRA while it is positive (with usually a larger magnitude) for 
NUGM:24. Given the relative weights of the A^ terms in Eq.5, it is thus quite apparent that 
the total contribution from the trilinear couplings to the evolution of B is quite similar in 
the two models. On the other hand, since the signs of are reversed in NUGM:24, the 
aforementioned cancellations in dB/dt would no longer be operative; rather, the different 
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contributions would enhance each other leading to a large AB. This is the very reason why 
we choose to concentrate on models like NUGM:24. We note in passing that although the 
RGE for B does not explicitly include the SU(3) gaugino mass, it implicitly depends on the 
latter via the contributions from trilinear couplings. 

We now discuss the dependence of B and Bo on tuq and the other parameters. Being 
obtained from the REWSB condition of Eq.4, B (and hence Bo) evidently depends on mo 
quite strongly. The structure of Eq.5 suggests that, to one-loop order, AB should not depend 
on m . However, a subsidiary dependence arises through the determination of the scale at 
which the minimizations of Higgs potential (i.e. REWSB) is to be performed. Canonically, 
this scale is determined by demanding that the contribution, to /i 2 , of the 1-loop correction 
terms of the effective potential be small. In our analysis this scale is approximately halfway 
between the lowest and highest mass of the spectra and, generally, is not very far from the 
average stop mass scale ^/m^m^ (see Ref. [34]). Since this scale does depend on m , it leads 
to a small dependence in AB as well by virtue of being a limit of integration for the RGEs. 

2.2 Incorporating CP violating phases: |i?o|/|-B| and phase natu- 
ralness measure $ 

Even on inclusion of phases for the A and B parameters, the RGEs formally remain the same 
as in Eqs.(5&6). The evolution of the phases can then be extracted by comparing the real 
and imaginary parts of the said equations. Clearly, unlike in the case of the real parts, the 
imaginary parts of the beta functions for A's and B do not depend on the gaugino masses 
and hence there is no cancellation between the different contributions. Furthermore, even a 
vanishing 9b can lead to a non-zero 9b provided Aq has a non-trivial phase. For example, 
in the small tan/3 limit, the explicit analytical solution gives 

\B\sin9B = |-Bo|sin#B Q 

|.B|cos#b = |-Bo|cos#b 

We examine now the interdependence between the phases, their evolution (also see Ref. 
[15]) and the phase sensitivity $ for different values of tan/? and other parameters both within 
mSUGRA as well as NUGM:24. As we have already mentioned, the EDM constraints limit 
9b to be tiny (^ 0.1, and typically much smaller). Now, if either of \Aq\ or <pA is small 



1 

3' 

1 

2 



1 - D )\A \ sm 4> Ao 

1 — _Do)|A)| cos(f>A — Cm i . 



(10) 
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(actually, if \A \ sin (j>A <C \B\ sin 9b), then 9b would be determined essentially by \B\, \B \ 
and 9b- In this case, <3a would be quite unconstrained. The dependence on tan/3 is crucial 
and is best understood by considering the two opposite limits, namely small and large values: 

• For a small tan (3 5 or so), sin 2(3 is large, and therefore \B\ is appreciably large (see 
Eq.4). Within mSUGRA, for not too large a value of \Aq\, the GUT scale value \Bq\ is 
then quite comparable to \B\. This can be understood by recognizing the cancellations 
between the various terms in Eq.8 that keeps C small and thereby keep B—B relatively 
small (courtesy Eq.7). Consequently, in such a scenario, 9b is not too different from 
9b- This remains true even for <pA = tt/2 which maximizes the EDM values [22]. 

On the contrary, the situation in NUGM:24 is quite different. Here, a larger difference 
between \B\ and \B Q \ is generated by the enhancement in C. Consequently, 9b becomes 
appreciably different from (and numerically larger than) 9b- 

• For a large value of tan/3, on the other hand, sin 2(3 is quite small. Thus, unless \y\ 
is extremely tiny (as happens, for example, in hyperbolic branch/focus point [34,35] 
scenarios), \B\ is constrained to be small and has only sub-dominant influence on the 
evolution of 9b- This, in turn, implies that the value of 9b becomes strongly correlated 
with that of 4>a ■ in other words, a high degree of fine-tuning in one will necessitate a 
similar degree of fine-tuning in the other. 

We now focus on the issue of phase sensitivity. As Eq.10 suggests, the range allowed to 
9 B (i.e. A9 B ) imposes rather strong limits in the 9b q -^>a plane. Adopting the measure of 
phase naturalness <J> (as espoused in Eq.l), one may estimate, from Eq.10, the amount of fine- 
tuning associated with the phase 9b - Now, as the RGEs suggest, the implicit dependence 
of $ on A occurs primarily through the dependence of B itself on A . Thus, to the leading 
order, one has an approximate relation of the form [15] 

$ ~ \B/B \ . (11) 

We would like to point out that although the above simplification (as also those of neglecting 
Db and y T ) is quite illustrative, we do not take recourse to it. Rather we solve the complete 
set of RGEs numerically and also compute $ numerically directly from its definition (Eq.l). 

Note that, as obtained from Eq.l and the first of Eqs.10, the measure $ actually involves 
a factor of cos#b in the denominator. This causes $ to be very large when 9b is close to 
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7r/2, as also a change of sign for $ when 9b crosses tt/2. We will see that this is indeed the 
case for NUGM:24 where 9b can easily cross tt/2 owing to a large degree of phase evolution. 
In the mSUGRA scenario, on the other hand, such a feature rarely appears. 

As we have already discussed, mSUGRA is associated with a relatively small degree of 
evolution in B, and hence \B\ ~ |-3o|- This leads to a low value of $ ~ 1 or, equivalently, 
to a high degree of fine-tuning in 9b ■ On the other hand, a non-universal gaugino mass 
scenario like NUGM:24 can provide us with a large evolution of \B\. This, of course, can 
generate either \B / B \ <C 1 or \B / B \ ^> 1. The parameter space corresponding to the 
latter case (which is typically satisfied better for smaller tan j3 zones) reduces fine-tuning in 
9b - We will see that the said reduction can be as large as a factor of 10 to 20 compared to 
mSUGRA. And finally, the very same large evolution of |JB| also implies that \B \ ~ could 
be a possibility within such scenarios. In NUGM:24 where the evolution of B is large, the 
above reduction of |So| toward zero is possible when \B\ is large i.e. when tan/3 is small. 
In mSUGRA too this is possible, but only to a limited degree, as the aforesaid evolution 
is smaller in extent. So \B\ needs to be closer to zero in order to have a tiny \B Q \. In this 
sense, a requirement of a smaller \B\ would then favor large values of tan/3 for mSUGRA. 
This we explore numerically in the next section. 

3 Results: Degree of 5-evolution and phase sensitivity 
for mSUGRA and NUGM:24 

We show our numerical results in two stages. To begin with, we examine the difference 
between the evolution of B in mSUGRA and the NUGM:24 scenarios in the absence of any 
phases. Building on the lessons drawn from this exercise, we investigate next the core issue 
at hand, namely the behavior of the phase naturalness measure $ in each of the scenarios 
and the differences therein. 

3.1 Results in the absence of CP violating phases 

Focusing first on mSUGRA, we begin with the value of B as determined, by the REWSB 
conditions, in terms of the other parameters of the model, viz, m , mi, A and tan/3. This 
study, coupled with that for the derived value at the GUT scale, B , would serve to indicate 
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the regions of the parameter space for which the phase sensitivity can be significantly reduced. 
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Figure 1: (a) TTie dependence of B and B$ on mi in mSUGRA with the other parameters 
fixed as shown. (b,c,d) The dependence of the ratio Bq/B on mi, m and A respectively, 
keeping the other parameters fixed. 



Fig. 1(a) shows the variation of B and correspondingly Bq with respect to mi. With an 
illustrative choice of parameters, viz. mo = 300 GeV, Aq = and \i > 0, we exhibit our 
results for tan 3 = 3 and 10. One finds that B, determined through the REWSB condition, is 
almost linear with mi. The dependence on tan/3, on the other hand, is quite nonlinear; but 
as already touched upon in the previous section, the REWSB condition implies that, for a 
given mi, B decreases with increase in tan 3. As for the evolution of B, we find that B ~ B 

2 

unless mi is quite large. This is reflective of the aforementioned cancellations between the 
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gaugino and trilinear terms of Eq.5 in mSUGRA. For our choice of Aq = 0, this is same as 
the cancellations between the terms of C of Eq.8. Once mi becomes large, the contributions 
from the gaugino part of Eq.5 dominates and the cancellations are no longer as effective. 
This causes Bq to supersede B as is shown in Fig. 1(a). 

The information regarding the evolution of B can also be parametrized in terms of the 
ratio Bq/B and this is displayed in Fig. 1(b) as a function of mi. This ratio is of particu- 
lar interest on account of its relatively straightforward relation with the phase naturalness 
measure $ (note that $ ~ \B / Bq\). As could have been guessed from Fig. 1(a) itself, the 
variation with mi is nearly monotonic. The shallow dip at small mi values is a consequence 

2 2 

of the variation in the degree of cancellation between contributions to dB / dt and is difficult 
to see analytically from the leading terms alone. For large mi, the ratio Bq/B is seen to 
increase with tan/3, while for small mi the behavior is opposite. This, within mSUGRA, 
indicates that a small value of mi, coupled with a large tan (3 seems to be best suited for 
achieving a low degree of fine-tuning in the phases. 

In Fig. 1(c), we display the dependence of the same ratio on m . While the behavior may 
seem intriguing at first, note that B depends on m only via the requirement of REWSB. 
As Fig. 1(a) has already shown us, for the reference value of mi = 300 GeV, B is typically 
somewhat smaller than B. Now, B grows smaller as mo decreases. Thus, for small mo and 
large tan/5, B can be very small and the aforesaid evolution implies that B would have 
been negative. On the other hand, for large m values, B is large and thus the relatively 
small evolution leaves the ratio Bq/B very close to unity. 

The dependence of Bq/B on the trilinear coupling parameter Aq is quite linear (Fig. 1(d)). 
This, again, can be deduced from Eq.7 where fixing tan /3, mi and m will give rise to a linear 
relation between Bq/B and Aq. Note that progressively larger values for tan/3 increases the 
importance of the trilinear term contributions to dB/dt, thereby increasing the slope of the 
curve. 

We now repeat the analysis for the case of NUGM:24 choosing Aq = as before. However, 
since the sign of the electroweak gaugino mass parameters are now reversed, the gaugino 
contribution to Eq.5 would now enhance the trilinear contribution instead of cancelling it. 
And since the sign inversion affects only the sub-dominant contributions to the evolution of 
At t b, the latter remain close to their mSUGRA values with the result that the total trilinear 
contribution to dB/dt suffers only a small relative change. The result is then a monotonic 
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decrease of B with an increase in mi , and hence, in an appreciably large amount of evolution 
(Fig.2(a)). 

A further consequence is that the ratio B /B too is monotonic in mi (Fig. 2(b)). The 
slope though decreases with mi, leading to a flat behavior for moderately large mi values. 
This can be understood by realizing that, apart from B being approximately linear in mi 
AB too is approximately linear especially for large mi. While the steep slope for small 
mi might seem intriguing given the almost linear behavior of both B and Bq in Fig. 2(a), 
it should be noted that B is very small for such mi and consequently any departure from 
linearity would be magnified in the ratio. That the slopes at small mi values grow with tan (3 
is understandable too, as for larger tan/5, the trilinear term contributions to AB assume 
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greater significance. 

The abrupt ending of the curves, especially for larger tan (3 values might seem curious. 
However, note that A T is appreciably larger in NUGM:24 than in mSUGRA (see Sec. 2.1). 
This leads to a rapid suppression of m^ , the mass of the lighter stau. While the latter also 
sees an enhancement on account of the SU(2) gaugino mass being significantly larger in 
NUGM:24 in comparison that within mSUGRA for an identical value of mi, this effect is 
sub-dominant. Consequently, for such parameter values, the lighter stau would have a mass 
smaller than the lightest of the neutralinos thereby becoming the lightest supersymmetric 
particle. Since this is phenomenologically unacceptable, such regions of the parameter space 
have to be discarded. Note though that the extent of the allowed parameter range in the 
mi-tan/3 plane does depend on the value of m . 

Fig. 2(c) displays B /B for different values of tan/3 as m Q is varied. As discussed before, 
B increases with increase of mo and diminishes with increasing tan /3. For most of the region 
(except when m is large and tan /3 quite small) the ratio can be large and negative because 
of a large degree of evolution of B in NUGM:24. For larger tan/3, B itself is much smaller. 
Hence a large evolution results into a large negative B . On the other hand, a larger value for 
mo pushes B higher and Bq would then be dragged down to a value near zero. Additionally, 
we like to clarify that the larger tan/3 curves really end near 2 TeV or so in Fig.2(c) because 
of the REWSB requirement. This is unlike the smaller tan/3 contours that span the entire 
m range displayed. 

As for the dependence on A (see Fig. 2(d)), the relationship is once again linear, as 
predicted by Eq.7, for either of the two models under discussion. 

3.2 Evolution of CP violating phases 

Having analyzed the simple case of 6b — 4>a = 0, we may now consider the effect of 
phases. To start with, we continue to maintain 6g — 0, but now consider <pA — tt/2, or, 
in other words, a maximal phase in the trilinear coupling. This choice maximizes the EDM 
values [22]. To study the generic features and compare with the results of Sec. 3.1, we first 
choose a relatively small value of \A \ (= 100 GeV). Thus, ^ft(Ai) and ^t(B) would not be very 
different from the analysis of Sec. 3.1 because of the absence of any phase in the gaugino parts 
of Eqns.(5 & 6) and the smallness of \Aq\. With this choice of inputs, the only contributions 
to dQ(B)/dt or d^s(Ai)/dt arise from ^(Aj) themselves, and hence there is no occasion for 
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cancellations/enhancements unlike in the case for the real parts. In addition, the effect of 
4>a on | .Bo | would be limited even for maximal 4>a unless \Aq\ is quite large. This is reflected 
by Figs. 3, wherein we display the variation of both \B\ and \B \ with mi for either model. 
The results are seen to be consistent with the no-phase cases of Fig. 1(a) and Fig. 2 (a). 




We now invoke a non-zero 9 b and analyze the resulting evolution of the same from 
the electroweak scale to the GUT-scale. In Figs. 4, we display this for both mSUGRA and 
NUGM:24, and in each case for two values of tan/3, namely 3 and 10. Again, for illustrative 
purposes, we choose, for the other relevant parameters, mo = 100 GeV, mi = 300 GeV and 
\Aq\ = 300 GeV with <pA — tt/2. Although the constraints from the EDM measurements 
restrict \&b\ to very small values (9(10~ 2 )), we display the functional dependence for a 
wider range of 9b- The apparent discontinuities for the NUGM:24 curves are not physical 
and have only been occasioned by the choice for the domain of 9b , namely [—it, 7r]. Clearly, 
the amount of phase evolution in NUGM:24 is seen to be higher than that in mSUGRA. 

Having established that the degree of fine tuning could, in principle, be smaller in the 
NUGM:24 case, we now perform a scan of the parameter space for both mSUGRA and 
NUGM:24 so as to quantify the extent of this reduction. In each case, we consider two 
different values of tan/3 (= 2, 10) while maintaining <pA = vr/2 so as to maximize the EDM 
values. Allowing mo, mi and \A \ to vary up to 1 TeV (with the lower end set in accordance 
with the current limits on super-particle masses), we show, in Figs. 5, the scatter plots in the 
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Figure 4: $b vs 9b for tan/5 = 3 and 10 with other parameters are as shown for mSUGRA 
and NUGM:24 scenarios. 9b (defined to lie in the range [— ir, ir]) is seen to be larger for 
NUGM.-24. 



$-mi plane. It is interesting to note that, for low to moderate values of tan/?, the measure 

2 

$ rarely becomes negative in the mSUGRA case, whereas in the non-universal scenario it is 
more evenly distributed. 

While |$| does tend to concentrate around zero (Fig.5(c)), note that, for small tan/3, 
the NUGM:24 case does have a significantly dense distribution up to |$| ~ 20 and values 
as large as |$| ~ 100 are also obtained, albeit with a reduced frequency. In contrast, the 
mSUGRA case barely registers a presence even for $ ~ 1.5 (Fig. 5 (a)). Thus, in going from 
mSUGRA to NUGM:24, the fine-tuning can be reduced by a factor as large as ~ 70. For 
the tan/3 = 10 case though, the improvement is much more moderate. As Fig. 5(b) shows, 
the mSUGRA scatter reaches up to $ ~ 3.5, whereas the non-universal scenario admits 
|$| ~ 10 (Fig.5(d)), or, in other words, a reduction of the maximal fine tuning by a factor 
of ~ 3. More important, though, is that the density of points at higher $ is much larger in 
the NUGM:24 case than for mSUGRA. In other words, it is far more likely to have a less 
fine-tuned point in the parameter space for NUGM:24. 

Concentrating on NUGM:24, we present, in Fig. 6, contour plots for $ in the mo — mi 

2 

plane for two different values of tan /3. Note that the limits on m and mi are 2 TeV, higher 
than what was chosen for Fig. 5. Once again, \A \ is fixed at 100 GeV with (f)^ = n/2. A 
comparison of the two plots clearly reinforces our earlier result that the fine-tuning is less 
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Figure 5: $ vs mi for mSUGRA and NUGM:24 for tan (3 = 2 and 10, when m and \A \ 
are scanned up to 1 TeV for <f)^ = tt/2. 



severe for low tan (3. Furthermore, the values of m and mi leading to a particular $ are 
highly correlated. Note that both signs for $ are possible. The region where $ changes 
sign is associated with a parameter point where 8b is ~ tt/2. To summarize, the results 
displayed in Fig. 5 and Fig. 6 show that it is indeed possible to obtain a surprisingly large 
amount of reduction of phase sensitivity even for relatively small sparticle masses. 

We now explore, in detail, the range of tan (3 that is associated with very low level of 
phase sensitivity or, in other words, a very large |$|. As has been argued earlier, \B\ itself 
strongly depends on tan/3. Moreover, AB, and thereby B too, has a nontrivial dependence 
on tan/3. Thus it is understandable that a very large |$| would indeed prominently highlight 
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|$| corresponds to a lesser degree of phase sensitivity. Switching of the sign of $ in some 
region of parameter space is associated with $b crossing it/2. Here \A \ = 100 GeV and 
<Pa = 7r/2. Smaller tan f3 cases have larger values for |$|. 



such a dependence. Rather than attempting a full, but very computing-intensive, scan over 
the entire parameter space, we choose to restrict ourselves to the subset of the parameter 
space that would naturally produce very large values for |$|, namely the region with small 
\B \ and small \A \. Hence we adopt a framework with given values for \B \ instead of 
tan/5. The requirement of REWSB determines tan/? once Bq, mo, mi and Aq are fixed. 
Note however, that the point Aq = Bq = would imply the absence of any SUSY CP phase 
at all scales. Thus, it is not surprising to obtain very large values of |$| in this scenario. 
However, in this part of our work the focus is simply to study, the effect of tan (3 on $ in 
detail, more importantly for large |$| values. To quantify our study of this issue, we choose 
small representative values viz. \Bq\ = 0.5 GeV and \Aq\ = 1 GeV, along with (pA = tt/2 so 
as to maximize the EDM contributions as before. In Figs. 7, we present various scatter plots 
for $ as m and mi are varied over a wide range (0 to 2 TeV). Note that the results of this 
analysis have a significant dependence on \A \. For example, increasing \A \ to 100 GeV may 
reduce $ by a factor of 10 to 20. As Fig. 7(a) shows, within mSUGRA, |$| could be as large 
as 100 while most of the points lie between 10 to 25. The situation is qualitatively different 
in NUGM:24 (Fig.7(b)) where |$| may go up to 1500 while typically ranging between 200 to 
600. Thus, NUGM:24 is much better able to accommodate low phase- sensitivity solutions 
than do the universal gaugino mass scenarios. 
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It is curious to note that, unlike what Fig. 5 suggested, $ could assume negative values 
within mSUGRA (see Fig. 7(a)). This prompts us present a scatter plot of $ against the 
derived quantity tan/3. As Fig. 7(c) shows, mSUGRA admits negative $ only for large tan/3. 
In fact, even for the positive branch, large values of |$| are typically concentrated in the large 
tan/3 (20 to 45) region. In contrast, for NUGM:24, $ assumes larger values typically for low 
tan/3 values (2 to 5). It should be remembered in this context that, within NUGM:24, the 
large tan/3 domain is significantly restricted from considerations of the LSP (see Sec. 3.1). 
That the favored range for tan/3 is different in the two scenarios is attributable to the 
interplay between the cancellations/enhancements in the RGE evolution of B on the one 
hand and the requirement of REWSB on the other. 

Finally, we comment on the case of \i < 0. It turns out that for this branch of /i and 
'pAo — tt/2, one has \B \ > \B\ for almost all the parameter space of NUGM:24. As a result 
one finds no advantage toward reducing the phase sensitivity. 

4 Conclusion 

As is well known, the experimental upper bounds on the electric dipole moments of the 
neutron and the electron impose strong constraints on any source of CP violation in su- 
persymmetric models, in particular on the weak scale phase parameters. For example, in 
the minimal supergravity model, 9b, the phase of the bilinear Higgs coupling parameter is 
constrained to be typically smaller than 0.01, with only some very limited regions (such as 
the focus point scenario) in the parameter space admitting slightly larger (^ 0.1) values. 
This, however, implies a severe fine-tuning condition for 9b , the value of the same phase 
parameter at the unification scale. In turn, <pA , the phase of the trilinear coupling param- 
eter is also severely fine-tuned. This has been a longstanding problem with mSUGRA-like 
scenarios. 

To quantify this problem, we define a phase naturalness measure <3> as the ratio of the 
spread of the phase 9b at the unification scale that is consonant with the spread 9b allowed, 
at the electroweak scale, by the electric dipole moment constraints A larger $ would imply 
a lower degree of phase sensitivity. One finds that, unless tan/3 is very large, $ may be 
approximated to B/B for much of the parameter space. 

In this analysis, we have demonstrated that models admitting a large RG evolution of the 
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Figure 7: a) and b): Scatter plots corresponding to large values of |<fr| /or mSUGRA and 
NUGM:24 cases. A smaller value for \B \ as well as smaller \A \ enhance |$| both in 
mSUGRA and NUGM:24- However, NUGM:24 is associated with much larger values for |$| 
as compared to mSUGRA. c) and d): Displays of the associated values of tan /3 vs $. tan/3 
is small (2 to 5) for NUGM:24, whereas the same for mSUGRA is large (20 to 45). 



bilinear Higgs coupling could be interesting in the context of a reduction in the fine-tuning 
of phases. In particular, we choose a supergravity-inspired scenario wherein non-universal 
gaugino masses arise from a gauge kinetic energy function f a p transforming as a particular 
non-singlet representation of 577(5) (NUGM:24 of Table 1). As in the mSUGRA (singlet 
fap) case, this representation, considered in isolation, introduces no additional phase for the 
gaugino masses. 
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Studying the nature of the evolution of B to understand the correspondence with phase- 
sensitivity, we identify the large cancellations in the RGE for B as being primarily responsible 
for the high degree of fine-tuning within mSUGRA. In the NUGM:24, on the other hand, the 
said cancellations are replaced by enhancements (on account of the reversal in the sign of the 
gaugino mass terms) and this translates into a reduction of the above-mentioned fine-tuning. 
In fact, $ can be significantly increased in NUGM:24 (by a factor of 10 to 20) with respect to 
comparable mSUGRA type of models. The said improvement is typically more pronounced 
for small tan (3 values. 

A particularly interesting result is the identification of extended regions in the NUGM:24 
parameter space which admit a low degree of phase-sensitivity even for relatively small super- 
particle masses. This feature is absent in mSUGRA as well as in most other models with 
high scale inputs for SUSY breaking. 

We further explored the dependence of our results, on tan /3, by specifically concentrating 
on the parameter space corresponding to very large $ (or very small phase sensitivity) so as 
to compare the two models. Naturally, this occurs close to vanishing A and B values. We 
adopt a scheme where B itself is given as an input parameter instead of tan/3, given the 
more direct relationship of Bq with $. Our analysis shows that, even here, the values of $ in 
NUGM:24 are typically larger by a factor of 10 to 20 in comparison to those in mSUGRA. 
And whereas mSUGRA generically requires large tan/5 (20 to 40) for |$| to be large, the 
NUGM:24 scenario prefers a smaller tan/3 (2 to 5) instead. 

Finally, while our analysis has focussed on SU (5) as the GUT gauge group, similar 
considerations hold for 50(10) as well. A suitable non-singlet representation resulting in a 
similar gaugino-mass pattern as in NUGM:24 would also produce such a reduction of phase 
sensitivity. 
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